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Abstract 

Dyonic classical solutions of SU(2) gluodynamics are discussed. 
Exact form of dyonic solutions in different gauges is presented and 
the nontrivial problem of composition of the dilute gas of dyons is 
settled. 

Classical interaction between (anti)dyons is considered both ana- 
lytically and numerically. Confinement in the dyonic gas is discussed 
in connection with the topological properties of individual dyon solu- 
tion. 

Fermionic zero modes of dyonic are displayed and the chiral sym- 
metry breaking in the dyonic gas is demonstrated. 
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1 Introduction 



The QCD vacuum is known to possess properties of confinement and chiral 
symmetry breaking (CSB). In absence of dynamical quarks (quenched ap- 
proximation) confinement is characterized by the area law of Wilson loop or 
zero average of Polyakov line, while CSB is connected to nonzero values of 
chiral quark condensate. It was found in lattice calculations [1] that both 
properties disappear at the same temperature T c , while a part of confining 
configurations survive for T > T c , ensuring area law for spacial Wilson loops 
[2] ("magnetic confinement" [3]). 

In addition, the QCD vacuum is characterized by the topological suscep- 
tibility x an d nonperturbative energy density or gluonic condensate [4] , both 
quantities imply (in terms of gas of topological charges, like instantons) a 
density of approximately 1 top charge per lfm 4 [5,6]. 

By now there is no model of QCD vacuum with properties of confinement 
and CSB, based directly on the QCD Lagrangian. 

The most elaborated model is the instanton gas or liquid model (IM) [6], 
which ensures CSB but lacks confinement [7]. Even so the instanton model 
shows rather realistic features for hadron correlators [8] demonstrating that 
CSB is already very important property for correlators. 

The same can be said about the Nambu-Jona-Lasinio model (NJL) [9] 
(not directly connected to QCD Lagrangian) where confinement is also ab- 
sent. Therefore in both IM and NJL hadrons can dissociate into quarks. 

Thus it is an urgent need to look for more realistic model of QCD vacuum 
which obeys both basic properties: CSB and confinement. The latter is 
associated widely with monopole-like degrees of freedom [10] , which may be of 
purely quantum or quasiclassical character. In the latter case one should look 
for classical solutions of Yang-Mills theory with monopole - like fields. These 
solutions are known for a long time [11]. Such solution can be obtained from 
the multiinstanton solution in the so-called 'tHooft's ansatz by the singular 
gauge tranformation [11]. In the general case of finite-action multiinstanton 
solutions the form of fields and generalized gauge transformation were found 
in [12]. 

The solutions have both color-electric and color-magnetic fields and we 
therefore shall call them dyons. 

The dilute dyonic gas has been suggested some time ago as a model of 
QCD vacuum and some simple estimates of Wilson loop has been done for 
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dyons of finite time extension [12], demonstrating nonzero string tension. 
Recently the interest for the dyonic solutions has revived. In particular 
lattice studies of a classical and quantum field of a dyon have been done and 
a qualitative quasiabelian picture of confinement due to dyons was suggested 
[13]. 

Meanwhile the CSB properties of dyonic gas has been studied [14]. It 
appears that each dyon has (infinitely many) zero fermionic modes [15] and 
therefore CSB may occur through the same or similar mechanism as in the 
IM, where each instanton has a zero fermionic mode and the gas of instan- 
tons create collectively chiral quark mass and quark condensate [16]. It was 
shown in [14] that this indeed happens, and the values of chiral mass and 
condensates depend on density of dyons in the gas. 

Recently in an interesting series of papers [17] the properties of the solv- 
able N = 2 SUSY, 4d model containing gluons, Majorana fermions and an 
adjoint Higgs field were studied. It was demonstrated [17] that there is a 
confining phase in the model with confinement driven by the condensation 
of dyons and magnetic monopoles. With that the dual-Meissner effect as 
confinement mechanism [10] obtains an independent support, and also the 
dyonic gas model gets an additional impetus. 

It is a purpose of this lecture to consider general properties of the dyon 
gas from point of view of confinement and CSB and to obtain estimates of 
the corresponding parameters. To this end we refine our ansatz of dyonic gas 
done 10 years ago [12] and try to obtain a self-consistent picture adjusting 
two parameters of the model: average size of dyon p and average distance 
between them, R. With some choice of p and R we show that one can reason- 
ably reproduce values of string tension a, quark condensate < qq >, gluonic 
condensate Gi = ^ < (F® u ) 2 >, topological susceptibility x- Moreover, one 
can calculate all field correlators and compare with existing Monte-Carlo 
data. 

As additional check of the model we study a possible scenario of temper- 
ature phase transition and discuss nonperturbative physics at T > T c . 

The lecture is organized as follows. In chapter 2 we define a single dyon 
solution, make a gas of dyons in chapter 3 and calculate Wilson loop average 
for a single dyon and dyon gas in chapter 4. We compare this result with the 
case of instanton and instantonic gas and demonstrate the reason why dyons 
confine while instantons do not. In computing string tension an important 
role is played by dyon-dyon and dyon-antidyon correlations and screening 
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length of dyon magnetic and electric charge, which we specifically study. 

In chapter 5 we calculate CSB for dyonic gas using resent results [14], 
and expressing effect in terms of p and R. 

In conclusion the deconfinement scenario for dyonic gas is discussed and 
some estimates for the dyonic gas are presented. 

2 Properties of dyonic solutions 

We remind the classical Yang-Mills solution in the so-called 'tHooft ansatz 

A a , = -- g % v d v lnW (2.1) 

where ff^ is 'tHooft symbol 

% v = e a^u, n,v=l,2,Z or Sou, /x = 4,or-6 aM , v = 4 (2.2) 

and 

Pi and Xi, i = 1, ... N are real; for finite pi = p one has the Harrington - 
Shepard solutions [11], but we are interested in the limit pi = p — > oo and Xi 
lying equidistantly along the straight line 

x k = f, kb, k = 0, ±1, ±2, ...±iVi, (2.4) 

In most cases the limit of Ni — > oo will be considered. 
For the choice p — > oo one has 

and due to (2.1) is in general periodic in time t. One can make a gauge 
transformation = A^f = U + (A^ + ^d^)U, with U = exp(i^6), such that 
A^ is time-independent for Ni — > oo [11,12]. In this limit one has 

= f(r)e iba n b (2.6) 
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Ma = ^{r)n a (2.7) 

with 

f(r) = -(l-^-) (2.8) 
gr sh^r 

1 2ir 

tp(r) = — ('jrcth'yr — 1), 7 = — (2.9) 

We shall call (2.6)-(2.7) the Rossi solution [11] or the dyon solution since it 
has both electric and magnetic fields; it is clear that Rossi solution belongs 
to the class of Prasad-Sommerfield solutions [11]. 

One can easily find color-electric and color-magnetic fields, E ka = —B ka 

B ka = 5 ak (-f - f/r) + n a n k (f> - f/r + gf) (2.10) 

The field (2.9) contains both long-range and short-range parts, indeed for 
large r one has f(r) ~ 1/gr and B ka ~ — S ^ L . To make separation more clear 
let us go from the hedgehog gauge(2.5-2.6) to the quasiabelian (or unitary) 
gauge, where the long-range part of B k is Abelian. This gauge rotation is 
given by the orthogonal matrix [18] 

ik = cos95 ik + (1 - cos&)viV k + sin9e ik iVi (2.11) 

where we have introduced unit vector e and defined 

cosO = en; Visind = —ei mn e m n n (2-12) 

One also has property 

OikUi = e k (2.13) 
Therefore the gauge transformed A'^, 

4„ = i Mb • ba - ^-0 lb d^0 lc e a bc (2.14) 

yields 

A' 4a = ^{r)e a (2.15) 

and 

B' ka = (-2f/r+gf 2 )n k e a +(-f-f/r)[cos65 ka +(l-cos6)v k v a -e k n a ] (2.16) 
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Now we notice that the long-range part - the first term on the r.h.s. of (2.16) 
- has a fixed color direction; by choosing e along the third axis, we have 



E> k = B' k (r^oo)~-±- 2 -n k -^ (2.17) 



of course this transformation from (2.10) to (2.16) not defined at 6 = ir. 
(Dirac string). We note also, that the long-range part of B k and E k can be 
written as a derivate 

B' = V<3> h short-range terms (2.18) 

2 

with 

$ = / (-2f/r' + gf)dr' (2.19) 
Jo 

The total action of dyon is proportional to its time extension 

S = \j d 3 fJ Q T dt{Bl k + El) = (2.20) 

It can be considered as a string of instantons of infinite radius [12], and for 
N centers in the string (the "N-string") one has 

S(N) = ^-(N-1) (2.21) 

The total number of parameters for the case when we allow the centers to 
move and have finite radii pi is equal to 5N — 1, (JV > 3) [19] (plus 3 overall 
color orientations). 

Thus the enthropy of the iV-string is smaller than that of N independent 
instantons - in the latter case the total number of parameters is 8iV due to 
independent color orientation of each instanton. This is the price one should 
pay for a new property of dyon: its coherent field is able to confine, as we 
shall demonstrate later. 

Now we turn to the quantum corrections around the dyon. There is a 
literature on the subject [20], but we shall use a simple method a la Polyakov 
to obtain the effective action of dyon due to quantum corrections around the 
dyonic classical solution. 
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One can write the effective action S e ff as [21] 



s en = \ J + ^ 2 ))f;m)f;a-^) (2.22) 



where the gluon plus ghost self-energy part n{q ) is 
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and one can introduce the renormalized charge 



1 1 11 1 Aq b n 2 , nni . 



g 2 (q) g 2 3 c 16tt 2 q 2 16tt 2 A 2 
For a long (anti)dyon AT ^> l) one can write 

± E ka = B ka = ((q )b ka (\q\) (2.25) 
with C(q ) = T si ^W , Introducing (2.24) in (2.21) one obtains 

s <» K T7WY (2 ' 26) 

where g 2 (q) is given in (2.24) and q is obtained from the integral / b\ a ' 

1 167T 2 



9 2 (<i) b 

2-i 4 

q « _L « (2.27) 

7T 

The form (2.23) is valid when g> A; for small q/A one should replace (2.23) 
by the expression which takes into account confining configurations in the 
vacuum [21,22] - dyonic gas in our case - and one has approximately 

1 ^ b _ ln l + <^ m 2 ^2na (2.28) 



g 2 (q) 16tt 2 A 2 
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3 Dyonic gas 



In this and following chapters the results of ref. [23] are largely used. Con- 
sider a system of several dyons and antidyons. We have several options for 
the composition of the system. 

(1) Dyonic g simple superposition ansatz (similar to the instanton 

gas ansatz [6]) 

N+ N- 

4.(*) = E^ (0 (*) + E4r (i) (*) (3.i) 
i=i i=i 

where N + and iV_ are numbers of dyons and antidyons respectively. Each 
individual vector potential Aft (with superscript + for dyons and - for an- 
tidyons) can be characterized by a Ad vector i?W and 0(4) unit vector uj^\ 
defining direction of the straight line passing through i?W and centers kb in 
(2.5), moreover there is an overall color orientation so that 

Af(x)=Qt(LA),(r,t)Q l (3.2) 

where A^(r,t) is the solution defined in chapter 2 and 

r=[(x- R®) 2 - (Or - R) ■ J l) } 1/2 , (3.3) 

t=( x -R)- W W , (3.4) 

L^ u (u) is the 0(4) rotation matrix, corresponding to uo. 

Thus the overall vector potential A^(x) in (3.1) depends on the set 
{Qj, R( l \ u®}, i — 1, N + + N_ and the total stochastic ensemble contains 

with all possible values of this set. 

For the case of zero temperature, QCD vacuum should be 0(4) invariant 
in the sense, that every observable K is an average over stochastic ensemble 
of the operator K(A) 

K =< K(A) >n,R,„ (3.5) 

where the weight of averaging is Poincare-invariant with respect to R, u and 
integration dfl is with the usual Haar measure. 

For nonzero temperature there appears a preferred direction in Ad. 

The relative simplicity of the ansatz (3.1) has a serious drawback -as 
also in the case of instanton gas - A^(x) is not a classical solution. As a 
consequence even two dyons has an interaction energy (action), and there 
are some additional divergencies, which we now discuss. 
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First of all one notices, that the action of the single dyon of fixed length 
L is finite together with quantum corrections and proportional to L. This is 
because B ka = ±E ka ~ \, r — > oo. For a pair of dyons situation is different. 
Let us consider the simplest case of dd or dd system with dyons in the same 
gauge (2.6-2.7), same u/ 1 ) = oo^ - two static dyons. 

A ia = e iba r lb H e iba r 2b (3.6) 

n r 2 

■4£ = ^k±^ 2o (3.7) 
r\ r 2 

and fi — r — is position of a dyon. 

For color magnetic and colorelectric field one has at large distances 

r W > 7 -i 

E ia = + EfJ + E% 2) , K ee Rf (3.9) 
where E^ and E 1 ^ is given in (2.10) and the interaction term E^ is 

E ^^l hS J^^±l) + r -2^± r JpL) (3.10) 
g T\r 2 r 2 r\ T\V 2 rfr 2 

It is clearly seen in (3.10) that the colorelectric contribution to the action 

S=lf(Ei + Bl)d 3 rdt (3.11) 

is diverging at large r for the system of two dyons, since E^ 2 \dd) ~ O(^), r — > 
oo; while for dd system one has 

E ia (dd) ~0(^), r^oo, i?= -i? (2) | (3.12) 

and the action is converging. From this it follows that it is impossible to 
have thermodynamic limit for N + ^ N_ (as well as for the Coulombic gas 
with nonzero net charge, here situation is similar). 

For the system of dyons and antidyons one can write for large distances 

N N h T 

i=l j^j=l z i/'j' j 
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E ka = f>£ + \ \-h a E ^(} + })+ (3.14) 

o r T 'i a T jk~n ( 

nrj rji-j J 

where 7$ = +7 and —7 for dyons and antidyons respectively. One can see 
that the only dangerous term in (3.14) is O(-) and it vanishes at large r if 
iV + = AL, so that Yh=i 1% — 0- Hence behaviour at large r for N + = AL is 
B ka ~ E ka ~0(l/r 2 ). 

Let us now calculate action of dd system, eq. (3.11) as a function of 
dastance Ri 2 between d and d. The leading term comes from (E\ 12) ) 2 and 
has the form 

S mt ~ X - J (E^) 2 d 3 fdt ~ T 7 2 it> 12 + 0(T 7 ) + ... (3.15) 



t 



One observes in (3.15) linear confinement between d and d, V int = ^f- 
I12. 

Even more striking is the behaviour of the d 2 d 2 system. One can compute 
total action of the system and Si n t = Stotai — Z)i=i 

Si n t depends on the configuration of the d 2 d 2 system when two d and two 
d are relatively close, while distance between d 2 and d 2 , R(d 2 — J 2 ), is large, 
one again obtains linear confinement 

S mt ~ Tj 2 R(d 2 - d 2 ) (3.16) 

For configuration of a "molecule" consisting of two atoms dd with large 
R(dd — dd) one gets instead 

r «~^ < 3 ' 17 > 

This behaviour reminds of a color Van-der-Waals potential between color- 
neutral objects. 

Thus a gas of static (3d) dyons and antidyons with N + = AL degener- 
ates at minimal action to the gas of dd atoms of size 7 -1 . It is intersting to 
compute fields and action of the dd system, when distance between dd, r(dd) 
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is not large. At r(dd) = the colorelectric field vanishes, Ei a = 0; while col- 
ormagnetic field B ia becomes short range, B ia (r(dd) = 0) ~ exp(-'yr), r — > 
oo. The total action is 

327r^ f°° dx 

S + T^— -tH(1 -xcthx) 2 + 2{l-x cos echx) 2 ] (3.18) 
g 2 b Jo sh 2 x 

Numerical estimate yields 

S(r(dd) = 0) ^ S jl, 1 (3.19) 

This should be compared with the action of each (anti)dyon, S = Soj. 

Thus one can see a strong attraction between d and d at small distances. 
The overal behaviour of the dd action numerically computed in [23] is shown 
in Fig.l. One can see both regimes (3.19) and (3.15). However, the lack 
of long-distance field, and zero topological charge of the dd atom makes it 
useless from the point of confinement, as we shall see in the next Section. 

The situation with the dyonic-antidyonic gas in 4d, i.e. when the (anti)dyon 
lines are have all directions, is different from the gas of static dyons. Indeed, 
let us take one dyon with the line along 4-th axis (i.e. static dyon) and an- 
tidyon with the line along the 3d axis. Assuming the superposition principle 
as in (3.6), one obtains for E\ a 

E£? ~ 7 W^nf (3.20) 

It is seen from (3.20) that the action (3.11) diverges for any distance between 
d and d, since (E^) 2 is nonzero at large r. 

Thus the 4d gas of d and d cannot exist in the ansatz (3.6-3.7) and one 
must choose another composition principle, which we consider next. 



(2) The 4c? model of dyonic gas 

We are using again the superposition ansatz (3.1-3.4). At this point one 
must specify in which gauge Aff(x) in (3.1) are summed up. 

We have found only one gauge (modulo global rotations) where the action 
of the superposition (3.1) is finite. This is actually the singular gauge of the 
original 'tHooft ansatz (2.1), 

4a = -^^hl^« 
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2r chr — cost 

yielding for a standard (not shifted and not rotated) solution the form 

m ,1 shr s „ sint 

A\a = e aik n k (- - cthr + 



1 r chr — cost chr — cost 

^ = „„(i_ dAr + _^_) (3.22) 



At large distances solutions (3.22) behave as 



A {i) ~ -, F {i) ~ — (3.23) 



and the same is true for the sum (3.1). 

Hence the total action of a gas of dyons of finite length is finite. The 
interaction energy of the dd system, V dd , and of the dd system, V dd in the 
gauge (3.22) depends on distance r and the relative time phase ip. At large 
r both V dd and V dd fall off as 1/r. This is shown in Figs.2 and 3. To ensure 
that dyonic gas could be a realistic model of the QCD vacuum one must 
investigate the following points: 

1) to check that the classical interaction between (anti) dyons is weak 
enough at large distances, so that the dilute gas approximation could be 
reasonably justified. 

2) to prove the existence of the thermodynamic limit for the dyonic en- 
semble (3.1), i.e. that the total action of the (big) volume V A is proportional 
to the volume, when it increases. 

One must also prove that the free energy of the dd calculated with quan- 
tum corrections has a minimum at a finite (and dilute) density. 



4 Confinement due to dyons 

We consider first the case of one dyon and calculate its contribution to the 
Wilson loop in two ways: 

i) first we use the long range part of the dyon field, appropriate for large 
Wilson loops as compared with the dyon radius ii) we show that Wilson loop 
can be rigorously computed through the function W of the 'tHooft ansatz 
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and find that the (magnetic) flux through the Wilson loop is n for the dyon 
and 2-7T for (multi)instanton; this explains why confinement is present for 
the first case and absent in the second in the simple picture of stochastic 
confinement [24]. 

Finally we evaluate the Wilson loop for the dd gas. 

i) Consider the circular Wilson loop as in Fig.4 and the dyon at the 
distance h above the plane of the Wilson loop (to be the (12) plane). We are 
using for a large loop of radius R the quasiabelian gauge form (2.17), which 
enables us to exploit the Stokes theorem 

W(C R ) = expig J F 12 da 12 = exp(-i^-ip) (4.1) 

where ip is the solid angle for the geometry of Fig.4, 

rR odo h 

One can see from (4.1)-(4.2) that for large R^> h (also R 3> b is necessary 
to use (2.17)) the color magnetic flux through the Wilson loop is 

flux(C) = \ = * (4.3) 

As we shall discuss later, this is the condition for confinement in the dyonic 
gas, in the picture of stochastic confinement [24]. 

ii) Consider now the general 'tHooft ansatz (2.1)-(2.3) and the same Wil- 
son loop, Fig. 4, where for simplicity we put h <C R, h — > 0. One has 

W(C R ) = exp ig J Adxi = expC^n^-) (4.4) 

where W is given in (2.3), and W r is derivate of W in r = \[x* at r = R. 

Two cases are possible; a) for the "periodic instanton of Harrington- 
Shepard [11], of the size p, when R^> p, one obtains 

W(C R )^exp(2t7TT 3 ^), (4.5) 

hence flux tends to zero (or 2n) and no confinement results. The same is true 
for one instanton [25], in the opposite case, p ^> R, again two possibilities 
appear, depending on the total lengh of the instanton chain L, L = Nib. 
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Indeed, for R ^> L, one has 



lim RW 

R^oo^ = -2 (4.6) 

and the flux is (— 2ir) - no confinement. 

Finally for the infinite instanton chain, L — > oo, 

W(R, x 4 = 0) = 1 SklR . - ^- (4.7) 

v ' ; 2R(ch~fR-l) 2R v ; 

and 

W(C fl ) = expHr 3 7r) = -1 (4.8) 

This coincides with the result (4.1) for h <C R. Thus the dyon cre- 
ates magnetic flux (— tt), while (multi) instanton creates magnetic flux (—2tt). 
Consider now the two-dimensional gas of d and d with 2d density j and with 
Poisson distribution w(n) 

w {n) = e- n ^f- (4.9) 

71/. 

The averaged Wilson loop can be written as 

< W(C R ) >= J2 e'" ^ 1 ^,^" = e ~ m = e ~° S ( 4 - 10 ) 

where the string tension a is : a — 2^. 

Thus the 2d Poisson gas of dyons (and/or antidyons yields confinement 
with string tension proportional to the (2d) average density of dyons. 

Next we consider the 3d gas dyons, which effectively means that in the 
total 4d ansatz (3.1) we keep for simplicity all dyons with roughly the same 
orientation = (0, 0, 0, 1). We again can use the quasiabelian gauge (2.17) 
and write for the dilute gas with N dyons inside volume V 3 

N J 

< W(C R ) >=< exp{-^-Y,^) >=< cos | > N > ( 4 - n ) 

1 i=i 1 

where 

ib f d 3 r ib Or) , „ _ „. 

< cos-^ >= / — < cos?-- 1 . (4.12) 



2 J V 3 2 
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Now for large volume V3, V3 R 3 , one has 



< W{C R ) >= exp(-aS) (4.13) 

where we have defined 

N ib 

a = — -In < cos— > (4.14) 
S 2 

For dyons distant from the Wilson loop one can write 

■kR 2 

ib(r) = — —cos6, (4.15) 

where 9 is the angle beween the direction to the dyon (vector r) and the 
perpendicular to the Wilson loop plane, R is the radius of the Wilson loop. 
Expanding in (4. 14) for small ib 2 <C 1 one has 

< >~ —;a « C— R (4.16) 

In (4.16) c is a numerical constant, iV-total number of dyons in the volume 
V 3 . Appearance of R in (4.16) actually violates the area law eq. (4.13). 
Indeed, a in (4.16) grows with the radius of the Wilson loop indefinitely; this 
situation may be called the superconfinement. 

One should note however, that the superconfinement occurs for the ideal 
gas of d and d, when one neglects completely correlations between d and 
d. In reality however for the tightly correlated pair dd the long range field 
disappears. Indeed, adding the vector potential (3.22) for the dyon and that 
of d, which obtains by changing the sign of A® and the second term of A\ z ^ 
in (3.22) one has 

— I s%ti hv 

A ia {dd) = 2e aik n k ( coth r H -) (4.17) 

r cos hr — cost 

A 4a (dd) = 

It is easy to calculate that the long range color magnetic and color electric 
fields for vector potentials (4.17) disappear, 

E ia (dd) = 0, B ia (dd) = 0(e~ r ) (4.18) 
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Hence for a given dyon, antidyons can partly screen its field and vice versa, 
and the phenomenon of Debye screening must take place. One can estimate 
the Debye screening mass for the Ad gas of dyons and antidyons [23] 

ml~£)W r D = ^- (4.19) 
V 3 m D 

so that the Debye radius is of the order of average distance between neighbor- 
ing dyons. This means that dyons which are farther away from the Wilson 
loop than m do not participate in the creation of string tension and one 
should replace in (4.16) R by r D for R 3> r D . Finally one obtains an esti- 
mate for Debye screened dd gas 

a ~ constm 2 D ~ A 2/3 (4.20) 

^3 

This behaviour of a ensures the area law for the Wilson loop, the supercon- 
finement due to the Debye screening transforms into the confinement. 

These estimates have been checked in [23] by numerical calculations of 
the Wilson loop. 

In Figs. 5,6 the contribution to the Wilson loop from a dyon or a tight 
dd pair is shown as a functions of position of d or dd. One can see that the 
dyon contribution is equal to n when dyon is inside the Wilson loop, and the 
dd pair contributes only when it is exactly on the Wilson contour. 



5 Chiral symmetry breaking in the dyonic 

gas 

We follow in this chapter the recent paper [14] . For the gas made of equal 

number of dyons iV + and antidyons AL, N + = iV_ = y in the big volume 

V4 we assume that the thermodynamic limit exists for the total action and 

other extensive quantities like the free energy, when N — > 00, V4 — > 00 and 

M- is fixed and finite. 
14 

We shall use for the dyonic gas the formalism similar to that exploited 
for the instanton gas by Diakonov and Petrov [5]. To study the CSB as man- 
ifested in the nonzero chiral quark mass and chiral condensate it is enough 
in case of instanton gas to consider only the case of one flavour, Nf — 1, 
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since the so-called consistency condition displaying CSB comes out the same 
also for Nf = 2, 3 [5,16]. Therefore we for simplicity confine ourselves in this 
section also to the case Nf — 1. 

The main driving mechanism for CSB is provided by the zero fermionic 
modes on the topological charge [6,5]. For the instanton case zero fermionic 
modes were found by 'tHooft [26], and later in [5] those have been used to 
demonstrate CSB in the dilute instantonic gas. 

In case of dyons fermionic modes have been found in [27]. For our pur- 
poses we consider two sets (They can be expressed one through another) of 
zero modes, one with continuous parameter (5 playing the role of quasimo- 
mentum. 

^W) = W i/i(d + id i a i )(W- 1 F^)U+ (5.1) 



where x = x 4 , 



gi/3n27T 



and U + is a constant spinor of positive chirality. 

Another set is labbeled by the integer n and is obtained from ip^ putting 
(3 — and keeping only one term in the sum over n. 

«,(*) = + ^m- \ 2+( "; +2m)2 ) (5.3) 

In what follows we shall use both sets. 

The main problem to be solved in this section is: given fermionic zero 
modes on each of dyons and antidyons; find the full quark Green's function 
for the dyonic gas with the superposition ansatz (3.1). 

To this end we make the same interpolating approximation for the one- 
dyon quark Green's function as in [5], i.e. in the exact spectral repre- 
sentation of % — 1, ...N 

s ^, y)=j: ^miM (5 . 4) 

„ X n — vm 



containing all modes n = 1,2, ...oo, we keep only zero modes and re- 
place the nonzero-mode contribution by the free Green's function, since they 
coincide at large n ~ \fp 1 - 
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Thus with So = (-iD(B) — im) 1 one has 

*«><*,») = *(*,»)+ E (5-5) 

, —im 

zero modes 

One can see that S {i) diver ges as m — > 0, we shall show however that the 
total Green's function is finite for m — > if iV + = iV_. 

Using (5.5) and (31) one derives the total Green's function to be 

s = s - 5>W(aO f-^->) . fc «8 )+ G/) (5-6) 

n,m 

where upper indices i, k run over all dyon numbers, 1 < i, k < N, while 
lower indices n, m run over all set of zero modes of the given dyon with the 
numbers i, k. We have also defined 

V}L = I u^{x)i{d - igB) U W(x)d 4 x (5.7) 

We keep here the field B to make the formalism gauge invariant; in estimates 
we systematically put B^ equal to zero. Note that and u m in (6.4) should 
have opposite chiralities, hence V lk refer to dyon-antidyon (dd) or opposite 
(dd) transitions, otherwise V tk is zero. 

One can introduce graphs as in [5] to describe each term in (5.6) as a 
propagation amplitude from a dyon i to a dyon k through scattering on 
many intermediate (anti)dyons centers, with scattering amplitude of each 
center (dyon) being 4^ and transition amplitude from center j (excited to 
the s-th level) to center / (excited to the r-th level) being VjJ. 

The lower indices are not the only new element in (5.5) as compared 
to the instanton gas model [5,16]. The zero modes u$ depend also on the 
Lorentz orientation of dyon, in addition to the color orientation fiW and 
position i?W Q f the dyon, see Eq. (3.2). 

uf(x) = n^u%\x - R®,u;®) (5.8) 

Our next task is to compute the matrix elements of ( — ^—^) lk fixing initial 

X im+V'nm 

and final states and averaging over all coordinates of intermediate dyons. To 
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this end we introduce as in [5] the amplitudes D lk m and P^ m for even and 
odd number of transitions V respectively 



In the definition (5.9) it is assumed that amplitudes of returns to the initial 
and final center % ar k are not included in D lk , P lk and should be added 
separately (which makes Eq.(5.9) not an equality, but rather a symbolic 
equation). This amplitude of the return to the center j we denote as 

A™ = D» n (RV\RV\nW,nV\uV\uW) (5.10) 

Since in integration over all intermediate coordinates (R( k \ Q( k \ lj^) is 
done, A mn does not depend on R^\Q^\u^ and is a constant matrix. 

Taking into account any number of returns to the same center j, brings 
about a matrix e mn , defined as: 



1 



£ 



mnn 

m 



-(1-imAU (5.11) 



With its help the equations, connecting P and D can be written as follows 

Pnm = ~h V ™h ~WI (5-12) 

im im 2v 4 J i 

D nm = -§-J *B®dn®duV±V«e,.P* (5.13) 

As a next step we separate out the dependence of P,D,V on lower indices 
and on Cl^\Cl^ k \ To this end we consider zero-mode solutions u$(x) in the 
form of (5.3) and make Fourier transform 

u W(p) = j u^ixy^cPx = e iP ° 2nn u^(p) (5.14) 

It is important that u^\p) does not depend on n altogether. Therefore with 
the help of (5.7) one has 

V^M^^-R^M*^) = J ^e^ R(i) - RU %ii m (p), (5.15) 
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v% m {p) = e- 2ni{p °- p o m h + ( P i )tt + ( i \-p)tt {j) u(p j ) (5.16) 
where p % = 3? Wi p, and 3?^ is 0(4) rotation transforming time unit vector into 

We introduce now "amputated" amplitudes d, /, w instead of D, P, V as 
follows 

D ™ = I -j^ eip(R(i) ~ R(k)) ~ i2 *^^ ( 5 - 17 ) 

and similarly for f{p\p k ); according to (5.16)) one has w{p\p k ) = —p 
Insertion of these definitions into Eqs.(5.12-5.13) yields 

d(p\p k ) = -^ w -w(pHp')d^ i 'f(p',P k ) (5.19) 
where we have introduced 

u(p) = e +ipo27rn u{p)e ns e~ ipo27rs u + {p) (5.20) 

n,s 

One can see in (5.18-5.19) that / and d do not depend on rotations in p\p k 
and the integration over did 1 there acts only on v{p>), so that with the defi- 
nition 

v(p) = J is( Pj )dLu j (5.21) 

one obtains 

iNpup 



d(p) = fff 2 (5.22) 

1 + (zvInJ P p P p 

and f(p) is expressed through d via (5.18). The definition (5.10) can be 
used now to obtain the selfconsistency relation, taking into account that at 
m — > 0, A ~ -kr and therefore one has 

^■mn^ns = n^ms (5.23) 

m z 
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as a result of insertion of (5.17) and (5.22) into (5.10) multiplied with e mn , 
one has 

_ 2V 4 N C r d'p MHp) 

n ° ~ ~i\r~ J (2^mh^t? (5 } 

where we have defined the average number of zero modes per dyon —no, 

n pa , b is the internal scale parameter of dyons, 
We also introduced the chiral mass M{p) 

where we used the fact that V is averaged over all directions and should be 
proportional to the unit matrix in Lorentz and color space. 

Eq.(5.24) goes over into the corresponding consistency relation for instan- 
tons [5.16] when hq = 1 and matrix e becomes a number, while u(p) is the 
Fourier transform of the 'tHooft's zero mode [26]. 

The solution d{p) (5.22) assumes the knowledge of the matrix e ns , while 
the consistency relation (5.24) imposes only one condition. Therefore the 
strategy of solution is as follows. From (5.17) one finds A mn = -D^„ = 
(f^e _27np °( Tra ~ n )'U 4 ~(p)d(p)u(p) through d(p). It clearly depends only on the 
modulus \m — n\. Then inverting the matrix A mn one finds e mn from (5.23). 
Finally from (5.20-5.21) one finds v(p) and inserts it into (5.22), defining d(p) 
The cycle is thus completed, and should be repeated till the convergence is 
achieved. 

One can also study another basis of zero modes, namely that of (5.1). 
In this case dependence on (5 can be also extracted, indeed 



u 



«(p) = Y j e~ l ^ n u ( i ) =J2e 27Tni{p °-P ) u (t \p) = 



£ S(P -po- fc)«W(p) = 5 [M u^(p) (5.26) 



k 



where we have introduced notation 5[/3 Po ] , implying that f3 is in the interval 
[0,1] and 5 - function should be moderated first, introducing finite number 
of centers N in the dyon (Y^n=-N /2) an< ^ considering limit N — > oo at the 
end. 
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In this way one obtains the same equations (5.17-5.19) for /, d if the new 
definitions are used, e.g. 

D%, = j P {2it)4 S^uHp^WpVHp^'po] (5-27) 
and where in (5.18-5.19) now V is defined as 

Hp) -»• Hp) = J d ^Hp) £ [ P0 M u+ (p) ( 5 - 28 ) 

and ^ 

£ \poM= J d @ f d/3'$lP,Po] £ W$\P',Po] ( 5 - 29 ) 

From (5.27) one deduces that App> = Spp>A(/3) and hence also Epp* is diagonal 
due to the relation 

i 

£ W A P>P" = Z^<V' ( 5 - 3 °) 

lit 



and is equal to 
with 



epp, = 6^e([3) = 5 p ^A-\[3) (5.31) 

III 



A(/3) = / J^Mo^HpHpHp) (5-32) 

The system (5.22), (5.24-5.29) is now complete. 

We now proceed to write down the quark propagator (5.6) in terms of 
functions d, f and finally in terms of the chiral mass M(jp) (5.25). 

Following the same procedure as in [5], one can rewrite (5.6) as 

s{p)= yMt + ( A T^A) ^ x 



ns i 



x J dn {i) du {i) (u^(p,u {i) )u^ + (p,u {i) )+ dyon <-> antidyon)- 

Af \ 2 r 

—J J dfi«dfi^)da;«^^\«(p,a;«)(m £ns )^K^^)^(P,^) 
N \ 2 r - — 
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+ (5.33) 

Submitting in (5.33) expression (5.17-5.19) and (5.22), (5.24) we finally 
obtain S(p) in the form 

This form justifies the meaning of M(p) as a chiral mass, i.e. an effective mass 
of quark due to CSB. It coincides with the form of S(p) for the instantonn 
gas [5], however the explicit expression for M(p) (5.25) differs. 

The most remarkable feature of (5.34) is the disappearance of the massless 
pole ^ from S (p) One should have in mind of course that the form (5.34) is 
gauge-noninvariant and obtained neglecting confinement. If one takes into 
account these effects, as in [27], the pole structure in (5.34) is supplemented 
by the area law due to the string between the given quark and an antiquark 
and the pole is never present in physical amplitudes. 

From (5.34) one can easily compute the chiral condensate: 

< qq >Mink.= -i < qq >Eud.= i < trS(x,x) >= 

■ f d 4 p f d 4 p M{p) 

= l l J2^ S{P) = - ANc J (27) V + M^) (5 - 35) 

It is nonzero thus confirming the phenomenon of CSB in the dilute dyonic 
gas. 



6 Conclusions and prospectives 

We have given arguments that the dilute dyonic gas provides confinement 
and CSB and therefore may be a good candidate for a realistic quasiclassical 
QCD vacuum. Additional numerical checks are necessary of the area law 
of the Wilson loop for the 0(4) invariant Ad dyonic gas, which are now in 
progress [23]. If confirmed, the dyonic gas ansatz can be used in the same 
program of detailed calculations as were done for the instanton gas [8]. In 
addition one can calculate field correlators and condencates to be used as 
input in OPE and the vacuum correlator method [3,22,27]. 
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Meanwhile in anticipation of exact numerical checks let us estimate roughly 
parameters of dyonic gas which could ensure realistic values of 1) gluonic con- 
densate 2) chiral condensate 3) string tension 4) topological susceptibility. 
To get 1) and 4) at realistic values one needs roughly density of one topo- 
logical charge per lfm A . This can be saturated by 3d density of dyons of 
1 dyon per lfm 3 and with b ~ lfm. Then the average size of dyon is 
= ^ 0.16/m and one expects from Eqs. (5.24), (5.25) and (5.35) to 
get a realistic (within a factor of 2-3) chiral condensate. Finally, with the 
given dyon density the string tension (4.20) will be of a reasonable order of 
magnitude, a ~ several units xfm~ 2 . Thus order of magnitude estimates 
show that realistic model of dyonic gas is feasible. 

As a last point in this lecture we discuss now a possible scenario of tem- 
perature phase transition in the dyonic gas vacuum. 

The confined phase is described by the gas of dyons - better to say, gas 
of dyonic lines which are oriented in all directions; at T = these directions 
(wjp, see Eq. (3.2)) are spread uninformly in the 0(4), but at T > the 
distribution of uo^ may be deformed. It is important, that d and d are not 
paired, i.e. the correlation length of a dd pair is of the order of average 
distance between D and d, or n~ 1//3 . 

The deconfined phase can be chosen in such a way, that all dyons with 
lines directed along axis 1,2,3 are paired, i.e. d and d form neutral dd atoms 
with average size of 7 _1 ; producing no long - range field. Therefore string 
tension in the planes (14), (24) and (34) vanishes - there is no confinement 
in the usual sense. 

However, dyons and antidyons with lines uj^} along the 4-th axis are not 
paired, the (id average distance is of the order of n -1 / 3 and the confinement in 
the spacial planes (1,2), (1,3) and (2,3) persists. One can distinguish between 
two phases. If in the confining phase one can write 

S^ f = j:(S u (d) + S u (d))=8S , S = ^ 

then assuming that the number of d and d does not change, in the deconfined 
phase one obtains 

gdeconf = £ + S4{d) + ^ 5 ^ 0) 

u=l 
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where we have put S u (dd) « 1.1. So in agreement with the estimate (3.19). 
Thus one can see that the gluon condensate changes by some 40% across 
the phase transition; this fact roughly agrees with the magnetic confinement 
model of ref. [3], yielding reasonable estimates of T c (see lecture "Hot non- 
perturbative QCD" by the same author). 

One can see that the dyonic gas model may explain the deconfinement 
transition in a sensible way, however exact numerical computations are nec- 
essary to elaborate the detailed picture. 

This work was financially supported by the Russian Fund for Fundamen- 
tal Research grant 95-02-05436. 
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Figure captions 



Fig. I. The interaction energy, V(r) = S int /T, for a dyon and an an- 
tidyon in the static gauge, Eqs. (3.6-3.7) at the distance r vs r/b From 
Ref.[23]. 

Fig. 2. The interaction energy V(r,(p) =< S int > /T, in units of Sq/T 
where Si nt is averaged over the time period b, for the dd system as a function 
of distance T^d/b = r and relative time phase tp. The logarithmic singularity 
at r = 0, ip = is cut off by hand. From ref. [23]. 

Fig. 3. The same as in Fig. 2 for the dd system. The absolute minimum 
of V dd ~ is at (p = ir, r = and is equal to —1.3S /T. 

Fig. 4. The Wilson loop of radius R in the (l,2)-plane and a dyon at the 
distance h above the plane. 

Fig. 5. Contribution to the Wilson loop from the dyon placed at distance 
h = zR above the plane of the loop and at distance r.R from the center of 
the loop vs z and r. From ref. [23]. 

Fig. 6. Contribution to the Wilson loop from the tight dd pair. Notations 
are the same as in Fig.4. From ref. [23]. 
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